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Abstract. Using the recently found hA amplitude, single and double inclusive jet 
production rates on the nucleus are studied in the perturbative QCD with N c — > oo. 
Jet multiplicities are found to grow with A and energy as A 2 / 9 s A where A is the BFKL 
intercept. Long range correlations are found to behave as A 4//9 s 2A and change from 
negative to positive at a certain high energy. 
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1 Introduction 

In the colour dipole approach of A.H.Mueller [1,2] with N c — > oo the interaction with 
a heavy nucleus is exactly described by the sum of fan diagrams constructed of BFKL 
pomerons, each of them splitting into two [3,4]. The equation for this sum [3-6] has 
recently been studied perturbatively in [7], by asymptotic estmates in [8] and finally 
solved numerically in [6]. The results indicate that the cross-section of a particle on 
the nucleus saturates at high energies to its geometrical value 2nR\ corresponding to 
the scattering on a black disk. In itself this is not a very surprising result: since the 
days of the old Regge-Gribov theory it has been known that the supercritical pomeron 
fan diagrams with a triple pomeron vertex lead to a constant cross-section [9]. The 
new elements in the recent developments is that it has been understood that in the 
high-colour limit no splitting of the pomeron in more than two occurs and that the 
fan diagrams seem to be unitary by themselves so that they indeed represent the full 
scattering amplitude on the nucleus. 

An interesting byproduct of the numerical solution in [6] is the found gluon density 
of the nucleus in the combined momentum-impact parameter space 

dxG(x,q,b) iV c , , m 
d^b =2^ h{y ' q ' b) W 

where the rapidity y is related to x by 

y = in- 

x 

Function h defined via the sum of BFKL fan diagrams below (Section 2) and found 
in [6] numerically is a soliton wave in y — In q space moving towards higher rapidities 
with a constant velocity and preserving its nearly Gaussian shape. A fit to numerical 
data gives 

h(y,q,b) = h e~ a( -^°^ b))2 (2) 
where £ = lng, h ~ 0.3 and a ~ 0.3 are universal and £ is nearly linear in y: 

£ = c{b) + A y, A = 2.3a (3) 

where we standardly define d = a s N c /ir. The term c{b) actually depends on a dimen- 
sionless parameter B: 

B = na 2 s AT(b)R 2 N 

where Rn is the nucleon radius and T{b) is the nucleus profile function. Numerical 
results for c(b) are well described by 

c(b) = -3.11 + (2/3) In B 

Actually the gluon density does not really show itself in the total scattering cross- 
section at high rapidities. As mentioned, the cross-section becomes purely geometrical 
and insenstive to any dynamics. One expects however that the gluon density will be 
directly felt in production processes, where it is responsible for hard collisions giving 
rise to observable jets. 
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In this paper we study simplest of the production processes: the single and double 
inclusive jet production in hA collisions. It has to be noted that these observables were 
also studied for the fan diagrams in the old Regge-Gribov theory. The most important 
conclusion which was drawn from this study is that the multiplicities are strongly 
damped, as compared to the naive eikonal result. For hA scattering one obtains a 
multiplicity independent of A (instead of ~ A 1 / 3 ). As we shall see this result does 
not hold for the BFKL fan diagrams: the jet production rate grows with A although 
somewhat slowlier than in the eikonal approach, approximately as A 2 / 9 . The energy 
dependence remains esentially the same as in older studies: the production rate grows 
with rapidity Y essentially as exp Ay, where A is the BFKL intercept. With a s ~ 0.2 
its value is close to 0.5, so that the multiplicities reach very high values at rapidities 
in the region 20 -r- 30. 



2 Single jet inclusive production 

In the BFKL framework, at fixed impact parameter b, a single scattering contribution 
to the forward ampitude for the interaction of the projectile particle with a nucleus 
has a form 

A 1 (Y,b)=tsg 4 AT(b) J d 2 ri d 2 r 2 p(r)G(Y, r x , r 2 )p N (r 2 ) = lis J dV^n)^, r x , b) 

(4) 

Here p and p^ are the colour densities of the projectile and the target nucleon re- 
spectively. Function G is the forward BFKL Green function [10] taken at an overall 
rapidity Y 

G(F ' ry) = 32^JL em " Loo [u 2 + (n- 1)74] p + (n + 1)74] W*"' (5) 



n=— oo 



where and <fi' are the azimuthal angles and 

u (v) =a(>(l) -Re^(l/2 + w/)) (6) 

are the BFKL levels. Due to the azimuthal symmetry of the projectile colour density 
one may retain only the term with zero orbital momenta n — in (5). The correspond- 
ing inclusive jet production rate 

is obtained from the imaginary part of (4) divided by s after the substitution 

G{Y, n, r 2 ) - J d 2 rG( Vl , r u r)V k (r)G(y, r, r 2 ) (8) 

where r's are relative distances between the gluons and Vk(r) is a vertex for the emission 

am a <_ , 

V k (r) = A A (9) 

local in r (a differential operator, the arrows shows the direction of its action). The 
rapidity of the produced jet relative to the projectile is y\ = Y — y. 
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As mentioned, in the limit N c — > oo the total scattering amplitude which includes 
multiple interactions in the target is given by the sum of fan diagrams shown in Fig 
1. To obtain the corresponding inclusive cross-section, one has to make the substiti- 
tion (8) in one of the pomerons. However due to the well-known AGK rules all such 
substitutions in the pomerons below the upper splitting point (Fig. 2a) cancel, since 
the other branch, counting from the upper point, can be both "cut" and "uncut" (i.e. 
represent both real and virtual processes). The only left contribution is the one with 
the substitution (8) in the uppermost pomeron (Fig. 2b). It is clearly seen from Fig. 
2b that the vertex V k is coupled from below to exactly the full amplitude <&(y,r, b) 
which multiplies the projectile density in Fig. 1 and substitutes the single pomeron 
exchange contribution $ 1 (Y,r 1 ,b) in (4). Thus one obtains the inclusive cross-section 
for jet production on the nucleus at fixed b as 

I(y, k, b) = 2 J d 2 r'd 2 rp(r')G(y u r', r)V k (r)$(y, r, b) (10) 

The total inclusive production rate as a function of y and k is obtained from (10) after 
the integration over b 

Function <p(y, r, b) = r, b) / (2irr 2 ) , in the momentum space, satisfies a nonlinear 
equation [6] 

?^A = -H<l,(y, q ,b)-^(y, q ,b), (11) 

where y = ay and H is the BFKL Hamiltonian for the so-called semi-amputated 
function [10]. As mentioned in the Introduction, this equation was solved numerically 
in [6]. 

The Laplace operator in (9) acting on $(y,r, b) gives 
V 2 Mv,r,b)=2nVyf ^ e ^(f,(y,q,b) = -2nV 2 J ^ e ^V^(y, q,b) 

= 2tt J ^ e *VV^(y, q, b) = 2tt%, r, b) (12) 

Note that in the integration by parts used to obtain the second expression one has to 
be careful, since function <f>(y, q, b) has a logarithmic singularity at q = of the form 
— \nq. As a result, action of the Laplacian operator in q generates a S 2 (q) term which 
in its turn leads to a constant in $ (actually unity). However this term is annihilated 
by the Laplacian operator in r, so that no extra term arises from the intergration by 
parts. 

Function h(y, q, b), which is the gluon density, up to a trivial factor (Eq. (1)) is just 
a Fourier transform of h(y, r, b): 

h(y,q,b) = q 2 V 2 q <i ) (y,q,b) (13) 

So we find the inclusive cross-section in the form 

I(V, k, b) = J d 2 r'd 2 rp(r')V 2 r G(yu r', r)e irk h(y, r, b) (14) 

We observe that the definition of h(y, q, b) as the gluon density introduced in [6] is 
supported by its role in the jet production on the nucleus. 
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Applying the Laplace operator to the Green function (5) (with only n = terms 
retained) we obtain 

r' r°° dup yiw ( u "> 

^ G (y^--^LwTijw {r ' /r) (15) 

We are going to study the production rate at y\ » 1, that is far from the projectile 
fragmentation region. Then we can use the asymptotics of (15) at high y\, which comes 
from v ~ 0. Presenting io[y) in the standard manner 



where 

we get for (15) 



u(u) = Ay - [3v 2 (16) 
A = a41n2, &nd(3 = aU((3) (17) 



V^y.r^e^ (18) 

With (18), integration over r' converts r' into the average transverse dimension of 
the projectile (r) P = R P . Doing the integration over r we obtain 



I(y, k, b) = -^Rpe^ Jjjjj-Hv, ^ b ) (19) 



mk , b) .f^ Hy ,r, b) .J^ (20, 



where 

d2r jkrui„. „ u\ _ f d2( i h (y,q, b 

2ir \k + q\ 

As mentioned in the Introduction, the gluon density h(y, q, b) has a sharp peak 
at q — qo(y,b) = exp£ (y,b) and is practically zero everywhere else. With growing y 
qo(y, b) grows exponentially, so that at high y we can safely neglect k in the denominator 
of (20) as compared to q. Then F becomes independent of k and given by 

^Lh(y,q,b) = j ^ dtqh(y,q,b) (21) 
where £ = In q. One easily finds that 

Jd£h(y,q,b) = l (22) 
Indeed Eq. (13) can be rewritten as 

h(y,q,b) = -^ ( j ) (y,q,b) (23) 

Putting (23) into (21) and taking into account that at £ — > +oo = const — £ and at 
£ — > +oo — > we obtain (22). So (21) in fact has a meaning of the average gluon 
momentum in the nucleus (q). Since h has a sharp maximun at £ = £,o(y,b) given by 
(3) one expects that F(y,b) will approximately be given just by q (y,b). Numerical 
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integration gives more accurate values which are well fitted by the formula similar to 
(3) but with slightly different parameters: 



e Aiy 



F(y, b) = qi (y, b) = 1.25 2 / 3 — — , A x = 2.38a (24) 

RN^/y 

With this approximate form of F we find the final expression for the inclusive rate 
production (in (GeV/c)~ 2 ) 



, tolM) = 1 . 2 «^ [ ^ r(6) ^ 1 ^-^_^_ (25) 

where 

e = A - Ai = 0.39a 

and is relatively small. As we observe, the inclusive jet production is maximal in 
nucleus fragmentation region and slowly falls as one moves to higher y. 

Since the total hA cross-section saturates at large Y an immediate consequence of 
(25) is that the multiplicity grows like exp AY, that is, as the cross-section generated 
by a single pomeron exchange. This conclusion is also true for the inclusive production 
from the sum of fan diagrams in the old Regge-Gribov theory with a local pomeron. 
However this theory predicts a much stronger y-dependence, corresponding to (25) 
with e = A. 

As to the A-dependence, a novel element in the BFKL fan diagrams is factor 
[TralAT(b)R%] 2 / 3 . For a nucleus with constant density inside a sphere of radius Ra = 
A 1 / ? 'Rq its integration over all b leads to a factor 

Since the total hA cross-section at large Y saturates to its black-disc limit 2nR\ the 
jet multiplicity results proportional to A 2 / 9 : 



= wJ d%I{yXh) = L2A2 '^ eAr -"im^) (27) 

So unlike the fan diagrams in the old Regge-Gribov theory, the BFKL model predicts 
rising of the multiplicities with A, although somewhat weaker than according to the 
eikonal model (oc A 1 ^). 

Although the gluon distribution in the nucleus results radically changed as com- 
pared to the single nucleon, in particular, strongly damped at small momenta, this 
seems to have no effect on the ^-dependence of the jet production rate at moderate 
k « qi(y,b): it remains proportional to 1/k 2 and growing as k — > 0. The inclusive 
cross-section integrated over all k remains infinite, as for a single BFKL pomeron ex- 
change. Thus discussing the integrated cross-section we have to cutoff the spectrum 
from below by some k m i n - The formal expression for this cross-section can be trivially 
found integrating (19) over all k > k m i n . A simple estimate of its magnitude can be 
made in the logarithmic approximation in the large parameter q\(y, b)/k m i n . Evidently 

f d 2 k f d 2 q f d 2 k 



{27r) 2 k 2 F (y, K b)6(k - k mm ) = J -JLh(y, q, b) J J^J^Oik ~ k mm ) 
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1 f d 2 q 1 

~ 2tt J 2^ h ( y ' q 'ty ln ( q / kmin ") = 7^:(? ln (?/^mm)> (28) 

So, up to factor 1/27T, integration over A; leads to the average value of q\n(q/k m i n ) in 
the nucleus. Approximating it by qi(y, b) ln(gi(y, b)/k min ), so that for fixed b one finds 

i(y, ») - / - *~> - ») >« (29) 

with qi(y,b) given by (24). 

At high y the logarithmic factor reduces to Aiy. Its effect, opposite to the expo- 
nential, is to flatten the distribution I(y) as compared to I(y,k). The A dependence 
evidently becomes somewhat stronger. Numerical integration of I(k, y, b) with function 
h found in [6] gives results which more or less agree with this approximate estimates. 
At aY = 6 for lead (A = 207) we obtain multiplicities which in the region 2 < ay < 5 
very slowly grow from 1.13 10 5 to 1.26 10 5 . The found A dependence can be fitted 
by ~ A ' 26 . The large magnitude of the multiplicity is due to the exponential factor 
exp AY = exp(aF4 In 2). 



3 Double inclusive cross-sections 

Passing to double inclusive cross-sections and using the AGK rules we find two contri- 
butions schematically shown in Figs. 3 a,b. 

A simpler contribution Fig. 3a corresponds to emission of both gluons from the 
upper pomeron. To obtain it we have to make two insertions (8) in the upper BFKL 
Green function. We get 



I a (yi,ki,y2,k 2 ,b) = 



(2n) 4 da 



Fig.3a 



-dyid 2 ki dy 2 dd 2 k 2 d 2 b i 

= 2 J d 2 r > d 2 r l d 2 r 2 p(r')G(Y-y u r\r l )V k M)G(y 1 -y 2 ,r 1 ,r 2 ^^ (30) 

where we assumed that Y >> y\ >> y 2 >> 1 in the lab. system. 

The only new element in this equation is the second BFKL Green function with 
Laplacians applied in both coordinates. Using (5) we find 

^^^^/^(sr <3i) 

Combining this with (15) for the first Green function (with the integration variable Ui) 
we find an integral over r±: 

( d 2 r x ^rf- 2 ^ 2 ^ = 7r(h/2) 2 ^-^- T{\ + 1^-11;) 

J 1 V 11 J i{v x - v - i0) T(l - iv x + iv) V ; 

As before we limit ourserlves with the leading contribution as both Y — y\ and y 1 — y 2 
become large, which correspond to small values of both v and v\. Correspondingly we 
put v — v\ — everywhere except in the exponents u(vi)(Y — y±) and u(v)(yi — y 2 ) 
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and in the denominator of (32). We find that only the real part of (32) contributes, 
which is equivalent to taking 

ir8(vi — v) 



i{yx — v — iO) 

Using (16) we finally get for the double inclusive cross-section (30) 



Ia(yi,h,y 2 ,k 2 ,b) = lQnR P ^e A ^^^— -F(y 2 ,k 2 ,b) (33) 

where function F(y,k,b) was defined in Section 2 (Eq.(20)). A remarkable feature 
of this cross-section is that it depends only on the rapidity of the slowlier jet y 2 and 
does not depend on yi altogether. In the approximation (24) the dependence on y 2 
results the same as for the single inclusive cross-section, namely exp(— ey 2 ), so that the 
cross-section slowly diminishes as y 2 increases. 

Using approximation (24), integrating over b and dividing by a^X we get the corre- 
sponding multiplicity distribution 

M , m , k2 ) = 1 .2^6^ 7 -| | e— (34) 

where 7 was defined in (26). 

The second contribution (Fig. 3b) corresponds to emission of both jets just after 
the first splitting of the upper pomeron. To write it we first take the contribution from 
all non-trivial fan diagrams to $ [6] 

a 2 N r Y r 3 
*(yy,b) = -^ I dy' J j[d 2 x i 8 2 (x 1 + x 2 + x 3 ) 



x\x\ 



G(Y - y', r', x 3 )$ (y', x 1: b)$(y', x 2 , b). (35) 



where we denoted Xj, % = 1,2,3 the transverse coordinates of the triple pomeron 
vertex. According to the AGK rules the corresponding contribution to the absorptive 
part with both lower branches cut is twice (35) with the opposite sign. To finally obtain 
the double inclusive cross-section we have to insert into both <3>'s on the right-hand side 
verteces 14 (r) as was done in (10). In this way we obtain the second part of the double 
inclusive cross-section as 



g 2 N c r r Y r 3 

h(yi, h, y 2 , k 2 , b) = J d 2 r ' p{r')d 2 r x d 2 r 2 J dy J ]J d 2 x i 8 2 (x 1 + x 2 + x 3 ) 



-\^-G(Y - y', r', x 3 )G(y' - y u x u n)^ (n)^(y u n, b)G(y - y 2 , x 2 , r 2 )V k2 $(y 2 , r 2 , b). 
%i% 2 

(36) 

We again assume Y » y x » y 2 » 1 so that the rapidity of the splitting point y' 
varies from y\ to Y . 
We easily find 

-!V^(^'.-3)~/^"(|) (37) 
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Using this and representation (15) for the other two BFKL Green functions we find an 
integral over Xi, i = 1, 2, 3 

J( Vi ) = [ 17 d 2 x i 8 2 (x 1 +x 2 + x 3 ) — - — x7 2ivi Xo 2iU2 x7 2iU3 (38) 
J t=i ' X1X2X3 

This integral diverges at large x^s for real z/'s. In fact, presenting the 8 function as an 
integral over momentum q we find 

J = 2vr f[ 2~ 2w < ^ V , 2 ~ m \ I ^(„ 1+ , 2+ ,3) (39) 



1=1 



T(l/2 + iui) J q 3 



which does not exist for real v due to the divergence at q — 0. To overcome this 
difficulty we have to shift the integration contours in z/'s into the lower half-plane, 
using the fact that integrands are analytic in the strip —i/2 < Imz/ < +i/2. We 
choose to do it in a symmetric manner putting 

Vi = -i/6 + P h i = 1, 2, 3 (40) 

With such a choice we get an extra factor q in the integral over q, which then gives a 
S function: 

(41) 

Now we pass to the integration over the rapidity of the splitting point. The integral 
is done trivially: 

/ dy'e y ^ 1+UJ2 -^ = (42) 

where we denote U\ = w(vi) etc. Multiplying by the rest exponential factors we find 
the overall rapidity factor as 

JY -y\)u>\ + {Y -y 2 )u)2 _ AY-y\)u3+{y\-y2)i>]2 

~ (43) 

Note that each of the two terms depends only on two of the three v variables. Therefore 
integrating over z/'s we can always take as independent variables just those which are 
present in the exponents in (43). 

Having this in mind we again use the fact that all rapidity differences in the ex- 
ponents are large and use the saddle point method. Athough the initial contours of 
integration in the independent z/'s lie along the lines Imz/ = —1/6, the saddle points 
remain at v — 0, that is at P = (1/6)1 This means that independent z/'s will be close to 
zero, as usual. The only price we have to pay for the shift in the contours of integration 
will be the value of the dependent v at the saddle point. 

Take the first term in (43), for which v\ and z/ 2 are independent variables and 
z/3 is the dependent one. At the saddle point we shall have v\ — is 2 — and z/3 = 
— (1/ 6)i — V\ — v 2 = — (l/2)i. Accordingly we put v\ — v 2 — and z/3 = — (l/2)i in all 
terms except in the exponent and the singular T-function: 

/ / s \ r(l + iv\ + iiso) 1 , , 

r (1/2 - iu 3 ) = Tim + ivt) = ; ^ r 44 
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Again only the real part contributes, so that (43) actually gives 

tx8(vi + i/ 2 ) 

Using (16) and integrating over v\ and v 2 we obtain a factor 



A 6 )jp(2Y-y 1 -y 2 ) ^ 

and an extra factor r' due to the fact that u 3 = —(1/2)? at the saddle point. 

For the second term in (43) the independent variables are v 2 and is 3 which are zero 
at the saddle point. The dependent v\ is equal to — (1/2)2 at the saddle point. The T 
functions in (41) will give 7r<5(z/ 2 + z/ 3 ), so that the integration over v 2 and u 3 will give 
a factor 

1 A(Y-y 2 ) I 7T 

A 6 ^(Y-y 2 ) (46) 

and we shall have an extra factor r\ due to v\ = — (l/2)i at the saddle point. 

Other factors in (36) are straightforward. Combining all of them we find the ex- 
pression for the double inclusive cross-section corresponding to Fig. 2b in the form 

/.<»,*.,«,, fa, ») = ~e ^-»-> /=^=ZF( tefa , r ) 



(r 2 >,F(!fc,*,,6) - Rpe-* y ~*^ 2Y Y V * h( Vl M,bj\ (47) 
where 

(r 2 ) P = / d 2 rr 2 p(r) (48) 



The second term is evidently exponentially small as compared to the first and we 
can safely neglect it, having in mind all approximation already made. So our final 
expression is 

h(yi,k 1 ,y2,h,b) = 
> ^ (2y _ y;) Ffa,, fa, fa. r) (49) 

Recalling our expression (19) for the single inclusive cross-section we can rewrite it as 



I b (y u h, y 2 , k 2 , b) = 16h2(r) ^ - 2Y- yi -y 2 Hyu K b)I{m > **' 6) (50) 

After the integration over b with the approximation (24) we get the corresponding 
multiplicity distribution: 



^ yi ,h,y 2 ,k 2 ) = 1.44^/ 9 -^-^ 7 2 -^e 2 ^- 



7T 



-e(2/4-2/2) / 

51n2 i? 2 , ' fc?^ V p yi y2(2Y- yi -y 2 ) 

(51) 

Evidently at high Y the part dominates due to the factor exp 2AF as compared 
with only exp AY in /i a 
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Correlations are determined by the ratio 

R{Y, Vi, W = —, rr 52 

H(y 1: ki)n{y 2 , k 2 ) 

(evidently it does not depend on k\ and k 2 ). Neglecting the part /i a we find from (27) 
and (51) 

R(Y , m , V2) > <J^ M^WIM (53 ) 
15 In 2 Rp \j it 2Y — y 1 — y 2 

If the two jets, forward and backward in the cm. system, are taken symmetric, that 
is, 

yi = -Y + y, y 2 = -Y-y (54) 



Eq. (53) becomes 

R(Y,y) 



1 (r 2 ) 



2\ 

P 



15 In 2 R 2 P \ 7T \ Y 



£r-t\ (55) 



With the rapidity gap between the jet fixed and growing energy the term Ay 2 /Y 
can be neglected, so that the right-hand side of (55) becomes independent of the gap 
2y and proportional to \fY . This means that eventually it inevitably becomes greater 
than unity, so that the correlations become positive. The concrete threshold for this 
depends on the values of a and ratio (r 2 )/Rp. With the Yukawa distribution of colour 
in the projectile the last ratio is equal 2. Taking a s = 0.2 we get for (55) 

~ 0.195V 7 ! 7 (56) 

So with these values correlations remain negative up to quite high energies correspond- 
ing to Y < 25. 

Note that in the framework of the old Regge-Gribov theory, due to the assumed 
relations between the parameters and smallness of the intercept A, it was argued 
that the ratio R for the part //& is identically equal to unity, so that correlations are 
entirely generated by the part ji a [9]. In the present approach this does not hold. Al- 
though the factors exponential in rapidities are indeed the same in //& and the product 
fi(yi, ki)fi(y 2 , k 2 ), the power factors are not, so that their final F-dependence is dif- 
ferent. Also the relation between the coefficients is not fixed but depends on a s . As 
a result the ratio R for //& is different from unity and depends on energy and a s , so 
that the correlations become mainly determined by //& and thus much stronger than 
assumed in earlier studies. 



4 Conclusions 

We have studied the inclusive jet production off the nucleus in the perturbative QCD 
with N c — > infinity, in which the total hA amplitude is exactly given by a sum of fan 
diagrams constructed of BFKL pomerons. We have used the numerical results obtained 
for this sum in our earlier paper [6]. Our main results are the following. 

Our formulas confirm interpretation of function h(y, q, b) defined by Eq. (13) as the 
gluon density in the nucleus. In [6] this was done by comparison with the perturbative 
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expression for the structure function, the validity of which for large values of the gluon 
density can be questioned. 

The inclusive jet production rate is found to rapidly grow with energy in the same 
manner as the cross-section generated by a single pomeron exchange. This result is an 
immediate consequence of the saturation of the total hA cross-section, and was also 
found in earlier studies of the fan diagrams with local pomerons. 

However in contrast to these older studies, our jet production rate is A dependent: 
jet multiplicities grow roughly as A 2 / 9 . So the A dependence is somewhat reduced as 
compared to the eikonal result A 1 / 3 , but not so strongly as thought before. This re- 
sult may have direct consequences for particle production in nucleus-nucleus collisions, 
where we may expect the multiplicities grow as A 10 / 9 for identical nuclei, less than in 
the eikonal approach (A 4 / 3 ) but greater than in the older fan diagram estimates (A 2 / 3 ). 

The double inclusive jet production is dominated by the contribution from Fig. 2b. 
Unlike older studies, this contribution is not cancelled in the correlation function and 
leads to correlations which start negative at lower energies and then change to positive 
at a certain (rather high) energy. Their magnitude corresponds to a double pomeron 
exchange and so grows very large at high energies. 

In conclusion we want to remark that since the gluon distribution in the nucleus 
is strongly shifted to higher momenta, the notorious difficulty of the diffusion towards 
very small momenta seems to be overcome. This problem always raised doubts of the 
validity of the perturbative treatment of hadronic amplitudes: even if the initial wave 
function is centered at high-momenta, its evolution in rapidity inevitably introduces 
low momentum contributions, for which the perturbative treatment is not valid. The 
large nucleus seems to damp these small momenta contributions, so that the use of the 
pertubation theory becomes justified. 
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6 Figure captions 

Fig. 1. The amplitude for hA scattering. Lines denote BFKL pomerons. 

Fig. 2. Inclusive jet production from the lower (a) and upper (b) pomerons. 

Fig. 3. Double inclusive jet production from the upper pomeron (a) and immediately 
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after the first spliting (b). 
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